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A 3-D numerical model is developed to represent the oscillating thermal- convection induced in a 
cylindrical cavity filled with Ga-doped germanium melt with and without the presence of an external 
magnetic field. The model is developed based on the penalty-finite element solution of the equations 
describing the transport of momentum, heat and solutal element as well as the electromagnetic field 
distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. 
Numerical simulations are conducted to study the convection and magnetic damping effects as a function 
of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied 
magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex 
recirculating convection pattern in three dimensions and that an applied magnetic field can be employed 
to suppress this deleterious convective flow and both magnitude and orientation of the applied field are 
important in magnetic damping of the g-jitter induced convective flows. 


1. Introduction 

Microgravity and magnetic damping are two mechanisms 
applied during the melt growth of semiconductor or metal 
single crystals to suppress buoyancy driven flows so as to 
improve macro and micro homogeneity of the crystals. 
As natural convection arises from gravity effects, 
microgravity offers a plausible solution to reduce the 
convective flow to achieve a desirable growth condition. 
However, flight experiments indicated that residual 
accelerations during space processing, or g-jitter, can 
cause considerable convection in the liquid pool, making 
it difficult to realize a diffusion controlled growth, as 
originally designed, when experiments were conducted in 
microgravity l *\ 

G-jitter is a phenomenon associated with microgravity 
environment and can have both steady state and transient 
effects on convective flows 4,5 *. Measurements by 
accelerameters aboard space vehicles reveal that g-jitter is 
inherent with space environment and that it is three 
dimensional and random in nature and has a wide band of 
frequency 6 *. Recently, the effect of g-jitter on fluid flow 
behavior in microgravity has been a subject of intensive 
research. Substantial work has been published to 
investigate convection flows in single crystal growth 
systems that are induced by both steady state and transient 
gravity perturbations under microgravity conditions. 
Numerical modeling results showed that the frequency, 
amplitude and spatial orientation of g-jitter are all 
important in affecting the convective flows 7 ' 14 *. The 
analyses further indicated that the most detrimental effect 
of g-jitter comes from the low frequency component with 
a large amplitude, especially when it is aligned 
perpendicularly with temperature gradient, which often is 


along the crystal growth direction. For crystal growth 
applications, the effects of the high frequency g-jitter 
components may be neglected. 

Magnetic damping effects on the motion of an electrically 
conducting fluid originate from the interaction of the flow 
field with an externally applied magnetic field. This 
interaction induces an opposing Lorentz force that results 
in a reduction of melt flow velocities. The use of 
magnetic damping has been widely exploited in the 
materials industries to obtain more homogeneous 
semiconductors and metal crystals under terrestrial 
conditions. While magnetic damping on turbulent and 
thermal fluctuations on earth has been well 
established 15,16 *; relatively little is known about the 
magnetic damping effect on the flow field generated by 
low frequency g-jitter in microgravity. 

Because magnetic damping stems from a different origin 
than gravity, it may be superimposed over microgravity to 
further suppress the deleterious effects associated with g- 
jitter. This idea has been explored recently 15-17 *. Both 
analytical and 2-D numerical analyses have been 
performed. The analytical solutions were derived for a 
configuration of a parallel plate thermal channel subject to 
both the oscillating gravity force and a transverse 
magnetic field and were used to serve a dual purpose of 
obtaining an asymptotic behavior of the system and of 
providing a check for numerical development. Numerical 
analyses were also conducted for a 2-D Bridgman- 
Stockbarger system for the growth of Gadoped 
germanium single crystals. These analyses showed that 
magnetic damping can be applied to reduce the g-jitter 
induced convective flows and the reduction of the 
velocity depends on the g-jitter frequency and applied 
magnetic field strength and direction. 
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This paper describes a numerical investigation of 
magnetic damping on g-jitter induced fluid flow and mass 
transfer in a 3-D Brigdman-Stockbarger system for Ga- 
doped germanium single crystal growth under g-jitter 
conditions. For this purpose, a transient 3-D numerical 
model is developed based on the Galerkin finite element 
solution of magnetohydrodynamic equations governing 
the fluid flow behavior in the presence of time dependent 
gravity perturbations and an applied magnetic field. It is 
important to realize that g-jitter is inherently three 
dimensional and time dependent and thus a 3-D 
dimensional model is more appropriate to fully resolve 
the g-jitter driven flows with and without the present of a 
magnetic field. 


2. Problem Statement 

The problem to be considered is schematically 
represented in Figure 1 below. In essence, we consider a 
simplified 3-D model for the Bridgmann-Stockbroker 
growth of Ge-doped silicon single crystals where the 
phase change due to solidification is neglected and as 
such the top and bottom surfaces are assumed flat. The 
melt is fed from the above into the molten zone. 
Solidification takes place at the lower side and the single 
crystal is withdrawn from the below. For the Bridgman- 
Stockbarger system, the lower side marked by h 2 is 
insulated and the upper side wall is fixed at the same 
temperature as the top surface or melt inlet. In 
microgravity environment, g-jitter varies in time and 
orientation to induce convection in the molten zone. The 
magnetic field is imposed to damp the convective flows. 
Shown also in the figure is the cylindrical coordinate 
system used for the numerical analyses. 


Fig. 1 Schematic representation of the 3-D model for 
magnetic damping studies. 

The physical phenomena occurring in the system 
described above are governed by the Maxwell equations 
for the electromagnetic field distribution, the Navier- 
Stokes equations for fluid flow and the scalar transport 
equations for both heat and alloy element distributions. 
For systems like this, the magnetic Reynolds number is 
small. Also, we are primarily concerned with the flows 
generated by low frequency g-jitter components. As such, 
the induced magnetic field may be well neglected in 
comparison with the imposed magnetic field and the time 
variation of the induced electric field is also small. These 
considerations allow us to reduce the Maxwell equations 
to one equation for the electric potential 15 “ 17) . With these 
taken into account, the mathematical equations governing 
the fluid flow, temperature and species distributions and 
electromagnetic phenomena in the melt system may be 
written in the nondimensional form as follows, 



V • u = 0 (1) 
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V 2 C> = V • (B x u)= B • (V x u) (5) 

The boundary conditions for the problem represent the 
physical constraints needed to solve the problem. For the 
present problem, the no slip conditions are applied at all 
the solid walls. Part of the side wall is thermally insulated 
and the other part is at a fixed temperature, the top and 
bottom walls are at the fixed temperatures and the electric 
current density is zero on the solid walls as no current can 
escape from the system, because the walls are electrically 
insulated as well. These conditions may be phrased in 
mathematical terms, viz., 
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In the above equations, u is the velocity, p the pressure, T 
the temperature, g* gravity vector, O the electric potential, 
B the applied magnetic field vector, t the time, C the 
concentration, n the outward normal, k the solutal 
partetion coefficient between the solid and liquid phase, 
and n the outward normal. Also the following scales have 
been used in nondimensionalization, L=d> U=a/L 2 , 
to=cPla y P = pc? I cP and O -aB, where p is the density, a 
the thermal diffusivity and B the applied magnetic field 
strength. The temperature is nondimensionalized as T = 
( &-T m )/(T h -T m ) where 0is the dimensional temperature, T h 
the melt inlet temperature and T m the melting point. The 
nondimensionalization process also resulted in a set of 
nondimensional system parameters, which are given 
below, 
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with g being the earth gravity, J3 the thermal expansion, v 
the kinematic viscosity, D the diffusion coefficient of the 
species, V m the growth velocity, cr m the electrical 
conductivity and p m the molecular viscosity. 


3. The finite element solution 

The above equations along with the boundary conditions 
are solved using the Galerkin finite element method. 
Applying the standard finite element procedures to the 
above equations and carrying out involved algebraic 
operations, one has the following set of matrix equations 
for the unknowns defined at the nodal points, 

M U+ K(U)U + -C P M .Cp r U 
e y 


+ C B T — Fu 

(20) 

M r T+K r (U)T = F r 

(21) 

M c C+K c (U)C=F c 

(22) 

K < j,<I> + C < t > U = F <p 

(23) 

where now U is a global vector containing all the nodal 
values of u, v and w. C, T and <£> are global nodal 


concentration, temperature and electric potential vectors 
respectively. M is the mass matrix, K(U)=N(U)+K^ is the 
stiffness matrix, where is the diffusion stiffness 
matrix, and N(U) is the advection stiffness matrix. C p is 
gradient matrix, C B the buoyancy matrix and Fu the 
"force" vector which incorporates the Lorentz force term 
and the boundary conditions. In the temperature, penalty 
continuum and electric potential equations, the subscripts 
T, P and O refer to the corresponding sets of coefficients. 
The above matrix equations are solved using the 
successive substitution method and the time derivatives 
are approximated using the implicit finite difference 
scheme. Both variable time steps with automatic error 
tracking and time step automation and fixed time steps 
may be applied. The majority of the results presented in 
the report are calculated using the fixed time steps. 





Fig.2 Finite element mesh used for numerical simulations. 


4. Results and Discussion 

The finite element model developed above enables the 
prediction of electromagnetic field, velocity, temperature 
and solutal distributions in the 3-D cavity under the 
influence of both g-jitter and an applied magnetic field. A 
selection of computed results is presented below. All the 
results were computed using 8-node brick elements and a 
typical mesh for the calculations is given in Figure 2. We 
are particularly interested in the effects of transient g-jitter 
and magnetic damping. Thus the results presented below 
represent the situation where a steady state thermal and 
fluid flow condition is established by a steady gravity force 
in microgravity and then g-jitter is introduced to perturb 
the fluid flow field. To investigate the magnetic damping 
effects, the magnetic field is switched on at the same time 
the g-jitter starts to act. As stated earlier, g-jitter is random 
and 3-dimesnional in nature. However, to obtain some 
fundamental insight into the g-jitter driven flow and 
magnetic damping effects, we consider a typical 
component of a synthesized g-jitter force and a simplified 
magnetic field arrangement. The physical properties for 
the simulations are given in an earlier publication 20 *. 



Fig.3fiflT snap shot (t=40 sec) of 3-D convectTve flow 
structure generated by a g-jitter component (amplitude = 1 
X 10' 3 go, go the earth gravity constant and frequency = 0.1 
Hz) oscillating in the x direction. The time is* measured 
from the point at which the time varying g-jitter force starts 
to act on the fluid. 


Convective flows in this system are bound to be complex 
because of the interaction of the incoming flows and the 
convection cells generated by g-jitter. Figure 3 plots such 
a typical fluid flow structure induced by a time varying g- 
jitter force 40 seconds after it becomes effective. The 
initial thermal and flow field is generated by a constant g- 


l€r force that has been discussed in detail by Pan et af^. 
In this case, g-jitter is assumed to act in the x-direction 
only but varying harmonically with time following a sine 
function. It is clearly seen from the figure that the 
convection pattern is complex and g-jitter force strongly 
perturbs otherwise a plug flow field from the top to bottom 
that would result only from the inlet inertia. 

For the system under consideration, the Pr number is 
small, suggesting that the g-jitter driven flows have little 
effect on the temperature distribution. This has been 
confirmed in this 3-D numerical study and as in the two 
dimensional case studied earlier 20 * the thermal transport in 
the system is by conduction only. This holds true also 
when an external magnetic field is imposed. 
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Fig. 4 Velocity field in the x-z plane at different time steps 
(viewed from the positive y-axis, same g-jitter as in Fig. 3): 
(a) t=5 sec, (b) t=10 sec, (c) t=35 sec, and (d)=40 sec. 



Fig. 7 Magnetically damped velocity field in the y-z plane 
corresponding to Fig. 5 but with an applied magnetic field 
in the x-direction B=0.22 Tesla. 


Fig. 5 Velocity field in the y-z plane at different time steps 
(viewed from the negative x-axis, same g-jitter as in Fig. 
3): (a) t=5 sec, (b) t=10 sec, (c) t=35 sec, and (d)=40 sec. 
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Fig. 6 Magnetically damped velocity field in the x-z plane 
corresponding to Fig. 4 but with an applied magnetic field 
in the x-direction B=0.22 Tesla. 


Figures 4 and 5 illustrate the fluid flow pattern in the x-z 
and y-z planes cut at different times after the time varying 
g-jitter effect comes into play. The three dimensional 
nature of the flow field is clear. In the y-z plane (or y=0 
plane), there exists a g-jitter component oscillating time 
harmonically along the x-axis. The oscillating gravity 
force produces an oscillating recirculating convective flow 
field in the plane. This picture is very similar to the 2-D 
case where a similar g-jitter gravity force is acting 20 *. 

In the y-z (or x=0) plane where no g-jitter gravity force is 
present, the situation is quite different in that no 
recirculating flows are induced and there exists almost no 
noticeable change in flow pattern (see Figure 5). Detailed 
analyses showed that the flow pattern in the other planes 
cut along the z-axis is somewhat between the above two, as 
expected. 
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Fig. 8 Magnetically damped velocity field in the x-z plane 
corresponding to Figure 4 but with an applied magnetic 
field in the y-direction B=0.22 Tesla. 


Let us now turn to the magnetic effects on the g-jitter 
driven flows. Figures 7 and 8 depict the convective flow 
field in the x-z and y-z planes with an applied magnetic 
field oriented in the x-direction and with a field strength of 
0.22 Tesla. Comparison of figures 4 and 6 clearly 
suggests that the applied magnetic has a strong effect on 


the g-jitter driven flows. The strong recirculation induced 
by g-jitter is substantially reduced and recirculation is 
suppressed in the region right above the lower side. In the 
y-z plane, the velocity field becomes virtually a plug field 
and the disturbance resulting from g-jitter is almost 
completely damped out, as appear in figures 5 and 7. 

From the theory of magnetic damping, the opposing 
Lorentz force depends on not only on the field strength but 
also its orientation. Figure 8 plots the velocity field in the 
x-z plane with a magnetic field oriented in the y-direction 
but with the same field strength as in Figure 6. By 
comparing figures 4, 6 and 8, it is clear that the magnetic 
field in the y-direction seems to produce a relatively 
smaller damping effect on the fluid flow in the x-z plane, 
as evidenced by the relatively strong convective 
recirculation persisting in the entire cavity. There is a 
reduction by the applied magnetic field but the extent to 
which reduction occurs is apparent less than by applying 
the magnetic field in the x-direction. However, in the y-z 
plane, the fluid velocity field is basically reduced by the 
magnetic field to a plug flow field, very similar to that 
shown in Figure 7. 


(figure will be insert here) 

Fig.9 Magnetically damped velocity field in the x-z plane 
corresponding to Figure 4 but with an applied magnetic 
field in the z-direction B=0.22 Tesla. 

By orienting the magnetic field in the z-direction while 
maintaining the same field strength, the fluid flow field in 
the x-z plane becomes different as shown in Figure 9. In 
this case, the recirculating flow induced by g-jitter is 
almost completely damped out and flow becomes 
approximately a plug flow. Analyses further indicated 
that in the y-z plane, the flow field is reduced also and a 
velocity field very similar to that in Figure 5 is obtained. 
These results seem to suggest that for this particular 
configuration, the maximum damping effect is achieved 
when the magnetic field is oriented in the z-direction. 


5. Concluding Remarks 

This paper has presented a 3-D numerical model for the 
study of transient g-jitter induced convective flows and 
magnetic damping effects in a simplified Bridgman- 
Stockbarger system for the melt growth of Ga-doped 
germanium single crystals. The model was developed 
based on the Galerkin finite element solution of tranient 
Navier-Stokes equations, thermal balance and solutal 
transport equations, and the simplified Maxwell 
equations. The results show that g-jitter induced flows 
are complex and three dimensional, evolve in time and 
oscillate in approximately the same frequency of the 
acting g-jitter force. For the system studied, a moderate 
magnetic field can achieve adequate flow reduction. The 
orientations of the magnetic field also play an important 
role in the reduction of the convective flows induced by 
g-jitter. 
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